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Abstract 

Motived by the necessity of explicit and reliable calculations, as a valid contribution to clarify 
the effectiveness and, possibly, the limits of the Tsallis thermostatistics, we formulate the Two- 
Time Green Functions Method in nonextensive quantum statistical mechanics within the optimal 
Lagrange multiplier framework, focusing on the basic ingredients of the related Spectral Density 
Method (SDM). Besides, to show how the SDM works, we have performed, to the lowest order of 
approximation, explicit calculations of the low-temperature properties for a quantum d-dimensional 
spin- 1/2 Heisenberg ferromagnet with long-range interactions decaying as l/r p (r is the distance 
between spins in the lattice). 
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I. INTRODUCTION 



Green's functions (GF's) are currently used in many-body physics and their power and 
success are widely recognized There exists at present a large variety of methods and 
techniques for calculation of the GF's both in classical and quantum thermostatistics 0, 0, 
0,0,0]. In particular, the related spectral density method (SDM), originally formulated by 
Kalashnikov and Fradkin j^J for quantum many-body systems, is a powerful nonperturbative 
tool which allows a direct study of the macroscopic properties of interacting quantum and 
classical many-body systems 0, 0, 0] also involving phase transitions. 

Recently, the increasing interest in Tsallis' nonextensive thermostatistics 0], has stim- 



ulated a lot of works 



10 



111 ] on the extension of the GF formalism also to this Tsallis' 



generalized framework of the statistical mechanics with the aim to provide new and effective 
methods for dealing with realistic nonextensive problems. Along this direction, the two-time 
GF technique and SDM have been formulated in classical nonextensive thermostatistics in 
two our papers |TJ| with application to the Heisenberg spin chain with short-range interac- 
tions. Here we wish to present the extension of the same formalism in quantum nonextensive 

n 

thermostatistics by using the optimal Lagrange multiplier (OLM) representation [12|, focus- 
ing on the spectral density (SD) and its spectral decomposition for their relevance in explicit 
calculations. Then, in order to show how the method works in the nonextensive context, 
we apply the extended SDM to a quantum <i-dimensional spin-1/2 Heisenberg ferromagnet 
with long-range interactions decaying as r~ p (r is the spin distance and p the decay expo- 
nent) and explore, to the lowest order of approximation, the nonextensivity effects on the 
low-temperature magnetic properties of the model. 

In Sec. II we summarize some basic ingredients of the nonextensive quantum thermo- 
statistics in the OLM representation for next utility. Sections III and IV are devoted to the 
extension of the GF formalism and of the SDM in the nonextensive context, respectively. 
In Sec. V the extended SDM is applied to the above mentioned spin model. Finally, some 
concluding remarks are drawn in Sec. VI. 
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II. THE OLM REPRESENTATION OF QUANTUM NONEXTENSIVE STATIS- 
TICAL MECHANICS 



The Tsallis' quantum thermostatistics |9| is essentially based on the so called g-entropy 
(with Ub — h — 1) 

1 - Trpl 

where p q is the generalized density operator satisfying the normalization condition Trp q = 1. 
In this framework, the generalization of the statistical average (referred as q- expectation value 
or q-mean value) for the observable O is given by 

Tr \p q a 6\ 

Here, q is a real parameter (called the nonextensivity index) which measures the degree of 
nonextensivity. 

Working in the canonical ensemble, the statistical operator p a can be determined adopt- 

nn 

ing the OLM constraint prescriptions in the extremization procedure of S q H2|| . Then, 
assuming the Hamiltonian H of the system to have a complete orthonormal set of eigenvec- 
tors {\n)} with eigenvalues {e n }, the normalized probability p n = (n| p q \n) associated with 
the n th eigenstate is given by 

p n = Z~ x [1 - (3(1 - q) (e n - U q )}^ , (3) 

where 

Z q = Y,[^-^-q){e n -U q )}^ ) (4) 

n 

and the g-internal energy U q = (H) q is given by 

U q = Z- 1 J2p g n en. (5) 

n 

Here (n\ p q \m) = p%8 nm are the matrix elements of the operator p q q in the {\n)} representation 
and Z q = Tr [p^] = YlnPn denotes the pseudo-partition function in the normalized OLM 
framework (Z\ is not merely the extensive partition function but contains an extra factor). 
In the previous equations, (3 = 1/T and T is the thermodynamic temperature. It is worth 
mentioning that, for q < 1, the formalism imposes a high-energy cutoff, i.e. p n = whenever 
the argument of the power function in Eqs. ©-© becomes negative 



III. TWO-TIME g- GREEN FUNCTIONS AND g-SPECTRAL DENSITY 



In the quantum nonextensive thermostatistics, we define the two-time retarded and ad- 
vanced GF's (g-GF's) for two arbitrary operators A and B as 

G i ;} B (t,t') = -ze u (t-t')([A(t),B(t%) q , (6) 

where v = r,a stands for "retarded" and and "advanced", respectively. Here, 9 a (t — t') = 
—9(t'—t) and 9 r (t—t') = 9(t—t f ), being 9(x) the step function; [..., ...j^ denotes a commutator 
(rj = — 1) or anticommutator (rj = +1); X(t) (with X = A,B) is the usual Heisenberg 
representation of operator X at time t. Of course, for g = 1, the conventional formalism 
is reproduced. It is worth nothing that, for general operators A and B, one can develop 
the g-GF's framework equivalently with commutators or anticommutators. However, for 
fermionic or bosonic operators it is of course convenient to use in Eq. ©77 = — 1 or 77 = +1, 
respectively. 



It is now relatively simple generalize most of the 



aasic properties of the conventional 
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two-time GF properties Q, |2| in nonextensive context 

Within equilibrium ensembles, the functions G q v \ B (t, t') depend on times only through 
the difference r = t — t' and hence one can introduce the Fourier transform 



"+00 

G';L M = / <1tG% (r) (7) 



Besides, in strict analogy with the GF extensive formalism 
density (g-SD) in the cu-representation is defined by 



/+oo 
dte^([A(r),B] v ) q , (8) 
-oo 



in terms of which it is easy to write the spectral representation of the associated g-GF's 



g&M= r ~. A rf K L,_ . ( 9) 



where (—l) u stands for +1 if v — r and —1 if v — a. These functions can be analytically 
continued in the complex cu-plane and combined to construct the single g-GF G q ab (^) = 
f +oc ^L A i AB ^ } f complex uj with a cut along the real axis. 

As a next step, one can immediately obtain the spectral decomposition of the g-SD. 
Indeed, by using Eqs. from the definition (JBJ) we have the exact representation: 



a \ Pn J 



A nm B mn 5 ( (10) 



with u mn = e m — e n . Then, the spectral decomposition for G q u \ B (u>) follows directly from 
the (^-representation 0. In particular, for G q ab we have 

GqAB n = i £ \p q n + wtt £^ (") 

y nm 



In strict analogy with the extensive case |4|, |y, |8|, Eq. (|TT|) suggests that the real poles u mn 
of G q AB{fjj) represent the exact energy spectrum of undamped excitations in the svstem. 

A comparison of the previous relations with the corresponding extensive ones jj, (| sug- 
gests that the Tsallis' statistics does not influence the meaning of GF singularities, but 
drastically modifies the structure of the spectral weights introducing a mixing of energy 
levels. 



IV. THE (/-SPECTRAL DENSITY METHOD (g-SDM) 

Now we have all the necessary ingredients to extend the SDM in the Tsallis' formalism. By 
successive derivatives of A q AB( T ) = ([A(r), B] n ) = f*™ duA q AB{u)e~ tu)T /2ix with respect 
to t and than setting r = 0, one obtains the infinite set of exact equations of motion or sum 
rules for A q AB{w) 

—u™k qAB {u) = ([U£A,B] ri ) q1 (m = 0,l,2,...), (12) 

-oo 271 - q 

where the quantity on the left-hand side of Eq. (J 12)) is called the m-moment of A q AB{&) and 
Ljj (m = 0, 1, 2, ...) acts as L° H A = A, L l H A = [A, H}_, L%A = [[A, H}_,H]_ and so on. 

The set of integral equations (|T2*j) represents a typical moment problem which should 
determine exactly the unknown g-SD. Of course, in practical calculations we must look for 
an approximated solution of A q AB (w) , which captures the essential physics of the system 
under study, truncating the set (fT2~J) at a given order. 

Suggested by the exact spectral decomposition (JTUJ), as a suitable approximation for 
A q AB{^)i one can assume a finite sum of properly weighted 5-functions (the so called polar 
ansatz) 

n 

K qAB {u) = 2^^(00 -J q %) } (13) 
k=i 

where n is a finite integer number and the unknown parameters A^b and ui^Jb are t° be 
determined self-consistently solving the first 2n moment equations in the set ()12J) . 
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In several situations of experimental interest, also the damping of excitations may be quite 
relevant so that the polar approximation (JT3j) is inadequate. In this case, to take properly 
into^ account the finite life-time of quasi-particles, one can assume, as in the extensive case 
the so called modified Gaussian ansatz 



A 



qAB 



(k) 

qAB~ 



exp 



Li) 



(k) 
qAB 



/r 



(fc) 

qAB 



k=l 



71 



r 



qAB 



(14) 



where T^ B represents the width of the peak at uj 
with frequency co^J B is given by rf$ B = 
It is worth noting that, the g-SDM is quite general and can be easily applied to different 
models by introducing the appropriate g-SD. 



.(*) 



oOqAB and the life-time of the excitations 



I\ a r> under the condition T^J B / 



(fc) 

qAB 
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<C 1. 



V. THE QUANTUM SPIN- 1/2 HEISENBERG FERROMAGNET WITH LONG- 
RANGE INTERACTIONS: A g-SDM APPROACH TO THE LOWEST ORDER OF 
APPROXIMATION 

Th e q uantum d-dimensional spin-1/2 Heisenberg ferromagnet with long-range interac- 
tions jl3( is described in the k-space by the Hamiltonian 

H = "2^ £ J W ( S * S -* + ^-0 ~ hS o- ( 15 ) 

k 

Here h is the external magnetic field, N is the number of sites {j} of a hypercubic lattice 
with unitary spacing; Sk, and J (k) are the Fourier trasforms of the spin operators Sj, 
= SJ ± iSj and the exchange interaction J; L j = J / — Vj\ p (p > 0, J > 0), respectively, 
and the k-sum is restricted over the first Brillouin zone (1BZ). 

For spin model (|15|). the g-expectation value of the magnetization per site is given by 

J=l k 

Due to the operatorial representation (|15 j) - (jlfij) , in the framework of the g-SDM it is conve- 
nient to introduce the g-SD 

A qk (u) = J°° dte- ( [S+ (r) , SZ k ] , (17) 
6 



which satisfy the infinite hierarchy of moment equations 

ho ° du 



. )7r ^A gk M = <[L™5 k + ,5: k ]_> (? , (to = 0,1, 2,...). (18) 



To the lowest order of approximation, we can assume for A g k(^) the one ^-function polar 
ansatz 

A qk (uj) = 2n\ q]5L 5(uj - u qk ), (19) 

where the unknown parameters A g k and o; g k can be determined by solving the first 2 moment 
equations of the set (|18j) . Then, working close to q — 1 and with the near saturation condition 
(S^S^y.) ~ (S£) (iS'f.k) = ^ 2m g^k,o [3], we obtain for uj q ^ and m q the set of self-consistent 
equations (with A 9 k = 2Nm q ) 

u q , = h + m Mmq + (20 ) 

1 1 v— v 2m 

m„ = > q - 5 , 21 

2 iV Z. [1 _ /5(1 _ g) ^ k] - T _ 1 ' ^ > 



where 



I ^ (k) = £(^£W =E iz^!L£. (22) 



It is worth noting that in Eq. (|2()jl . which determines the g-excitation spectrum u; g k, the 
contribution h + Jfi p (k)m (? is formally identical to the known Tyablikov dispersion relation 
|l|, |2j for the corresponding extensive problem. 

The solution of the self-consistent problem ()20|) - (f2*2*|) is rather complicate and one must 
consider asymptotic regimes for obtaining explicit analytical results. For instance, in the 
low-temperature ferromagnetic phase, we can resort to the Tyablikov-like approximation 

uj qk « 4 T) = h + - Jfip(k), (23) 

as a zero order in Eqs. (|2U|) - (|21|) . So, in the thermodynamic limit iV — > oo, Eq. (|21|) yelds 
a cumbersome expression for the g-magnetization m q (j3, h) which, for q = 1, reproduces the 
extensive counterpart and can be used to obtain explicit expansion as q — > 1. 

Interesting representation for m q ((3, h) and the g-susceptibility x 9 (/5, /i) = dm q (fl, h)/dh 
in the nearly saturation regime can be obtained, for d < p < 2d, under condition 
(3h(q — 1) > 1. Assuming the low- A; expansion fip(k) « A d (p)F- d [ij, where A d (p) = 
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7T d d d -P \T(p)]~ d / sin[7r(p — d)/2] and T(z) the gamma function, we find indeed 



m q ((3,h) 



Xq (P, h) 



K d A d 



2 d(l + Ph(q-l))* 

K d A d (3 

{p - d) (I + [3h(q - l)) 2 -^ 



Fi 



d 



p . JA d (p)P(q - 1)AP 



p-d'q-Vp-d' 2(1+/?%-!)) 



q-1 



(24) 



1 + 



JA d (p)f3(q-l)AP-<i 



g-l 



+ 



i-(2-^; 



d 



2(l+f3h(q-l)) 

q p JA d {p)(3{q -l)AP- d 



(25) 



K p-d' q-V p-d' 2(1 +/3h(q- 1)) 

Here, 2-F1 («, c; z) is the hypergeometric function, K d = 2 1 ~ d ix~ d l 2 /T(d/2) and A is a wave 
vector cut-off related to the 1BZ of the spin lattice. 

If, additionally, we assume (3h{q - 1) > 1 as T = — > 0, Eqs. (|2^|) and (|23j) yeld 
m g (/3, /i) ~ 1/2 — ^(/j)!^ and x<? (A ^) — B^ityT^ , where the explicit expressions of 
•A? (ft) and B< 9) (/i) (with /i ^ 0) are inessential for our purposes. 

The intrinsically nonextensive region (p < d), which is not of primary interest in this 
short contribution, requires a more delicate analysis which will be the subject of a future 
work. 



VI. CONCLUDING REMARKS 



In this short note we have extended, in nonextensive quantum thermostatistics, the two 



;ime GF formalism and the related SDM already developed for quantum 



Q3D 



and classical 



extensive systems. This offers the possibility to explore, at least in principle, the 
properties of realistic systems by using the big amount of experiences acquired in extensive 
problems. In case of the Heisenberg model (|15J). the polar ansatz (|19|) yields reasonable 
results for the excitation spectrum and the relevant thermodynamic g-quantities in the low 
temperature regime. For describing other thermodynamic regimes in a wider range of tem- 
peratures, one could to adopt a new set of approximations involving additional decoupling 
procedures and higher order moment equations, consistently with the spirit of the SDM 

mm. 
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